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Abstract 

The studies of fluctuations of the one-dimensional Kardar-Parisi-Zhang univer- 
sality class using the techniques from random matrix theory are reviewed from the 
point of view of the asymmetric simple exclusion process. We explain the basics of 
random matrix techniques, the connections to the polynuclear growth models and a 
method using the Green's function. 

1 Introduction 

The one-dimensional asymmetric simple exclusion process (ASEP) [1,2] is a stochastic 
process of many particles which perform asymmetric random walks with hardcore exclusion 
interaction. The ASEP is a simple model of statistical mechanical systems driven away 
from equilibrium [3,4]. In spite of its simplicity, it shows a rich variety of phenomena and 
has attracted much attention. The stationary properties at a single time are already of 
great physical interest. A lot of insights have been obtained about the boundary induced 
phase transitions, the long range correlations and so on using the powerful matrix product 
method [5]. 

As for time dependent properties of the ASEP, the average behaviors are rather well de- 
scribed by the mean-field approximation. We can also understand some part of fluctuations 
from the fact that the ASEP belongs to the Kardar-Parisi-Zhang(KPZ) universality class 
of surface growth models [6-8]. In particular the critical exponents are understood by the 
dynamical renormalization method and can be confirmed by some exact solutions [9,10]. 
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In the past decade we have seen a substantial progress of understanding of this scahng 
behaviors of the ID KPZ systems. In the language of the ASEP, a different approach based 
on a connection to the random matrix theory has been developed for the infinite system. 
It has allowed us to compute several time dependent quantities quite explicitly. Several 
important results and insights which have been obtained so far are (l)The fluctuations of 
the integrated current for several initial conditions have been computed and are shown to 
be equivalent to those of the largest eigenvalue of Gaussian ensembles. (2) The stationary 
two point correlation function has been computed explicitly. (3) The effects of boundary 
and initial conditions have been made quite explicit for certain special cases. 

A problem is that these results are rather involved and often described in the language 
of random matrix theory. Moreover some results are stated from the point of view of other 
related models and the connection to the ASEP is not obvious. Hence for those who are not 
very familiar with random matrix theory and these models, the developments might look 
not easily accessible. The purpose of this article is to explain the basics of random matrix 
techniques and the relationship to a few of the related models. Some part of this paper 
will overlap with the paper by Prahofer and Spohn [11] which contains a lot of information 
on the same topic. The present article aims at giving some background knowledge for [12] 
and at the same time a collection of new results obtained after [12]. See also recent review 
papers on related topics [13-16]. 

The article is organized as follows. In the next section, we deflne the model. We mainly 
study the discrete time TASEP with parallel update scheme. The continuous time version 
is also considered. In section 3, we map the problem of current to a combinatorial one and 
get an expression in the form of multiple integrals. In section 4, the basics of the random 
matrix techniques are explained. The multi-point fluctuations are studied in section 5 
by introducing the discrete polynuclear growth model and its multi-layer version. Some 
generalizations and variants are explained in section 6. In section 7, we explain another 
approach based on the Green's function. The flnal section is devoted to the concluding 
remarks. 

2 Models 

In this article we mainly consider a discrete time TASEP (dTASEP) on the inflnite lattice, 
Z. There are several versions of the update schemes. Here we employ the parallel update 
for which the model is deflned as follows. Each site of the lattice can be occupied by a 
particle or is empty (Fig. 1). Suppose that at each time step each particle tosses its own 
coin which shows heads with probability 1 — q. If it turns out a head the particle tries 
to hop to the right neighboring site. If the target site is occupied, the hopping does not 
occur due to the exclusion interaction among particles. If it turns out a tail, on the other 
hand, the particle remains at the same site. This is the TASEP with parallel update. The 
parameter q here should not be confused with the asymmetry parameter in the partially 
ASEP which is often referred to as q in the study of ASEP. 

The above deflnition describes how the conflguration of the lattice changes as time goes 
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on. It is also possible to define the same process by focusing on when each particle makes 
jumps to the right. Notice that the distribution of the waiting time of each hopping after it 
becomes possible is geometrically distributed with parameter q. Suppose that each particle 
has its own clock instead of the coin. It rings after geometrically distributed time at which 
the particle tries to hop to the right neighboring site. If the target site is occupied, the 
hopping does not occur. This description of the process defines the same TASEP with 
parallel update as above. 

The stationary measure with constant density p is rather simple. It is determined by 
the condition that the probability P{rij,rij+i) that arbitrary consecutive two sites j,j + 1 
are occupied {rjj = 1) or empty {rjj = 0) is given by [17] 

P(00) = l-p- J/(l-g), (2.1) 

P(01)=P(10) = J/(l-g), (2.2) 

P{n)=p-J/{l-q), (2.3) 

where J is the average current, 

J=i(l- Vl-4(l-g)p(l-p)). (2.4) 

Wc remark that in the stationary state the waiting time of each particle before making a 
hop is geometrically distributed with parameter 1 — J/p. 

In the study of nonequilibrium systems, one of the most fundamental quantity is the 
current. Let N[t) be the number of particles which crossed the bond between sites and 
1 from time through t. The average behavior is easily seen to be N{t) ~ Jt, or more 
precisely, 

lim ^ - J. (2.5) 

t-*oo t 

Hence the next step is to understand the fluctuation of this quantity. From the KPZ theory 
the exponent of the fluctuation is expected to scale as N{t) — Jt 0{t^^'^). The random 
matrix techniques allow us to compute even the scaling function. Note that the stationary 
measure is not enough to get such information. 

There is an important limiting case; g — > 1. Since the hopping probability 1 — q goes 
to zero, one has to rescale the time at the same time to take a meaningful limit. Let us 
consider the limit g — > l,t — oo with t = (1 — q)t fixed. In this limit the process reduces 
to the continuous time TASEP (cTASEP) which is the most well studied in the study of 
ASEP [1,2]. In cTASEP during infinitesimal time duration dt, each particle tries to hop to 
the right neighboring site with probability dt under the same exclusion principle as for the 
discrete case. It is again possible to define the model using the waiting time for hoppings. 
The only difference as compared to the dTASEP is that each clock now rings after an 
exponential time with parameter one. A remark is in order. In the following, when the 
cTASEP is discussed, the time t should be understood as the rescaled time i. 

The stationary measure of the cTASEP with density p is BernouUi with parameter p, 
i.e., each site is independent and is occupied with probability p or is empty with probability 
1 — p. The average current is known to be J = p(l — p). 
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3 Combinatorial Approach 



In this section we explain a result for the fluctuation of N{t) following [18]. The initial 
condition is taken to be the step one, in which all sites x < are occupied and all sites 
X > are empty at time (Fig. 2). In this case one can label the particles from the right 
so that the particle j {j = 1, 2, . . .) starts from site 1 — j. Several results for other initial 
and boundary conditions will be explained in sections 6 and 7. We denote the measure 
corresponding to the ASEP dynamics by P. 

The first step in this approach is to introduce the "waiting time" table, or the matrix 
w = {wij}, in which the {i,j) element Wij is the waiting time of the jth particle before 
making the ith hop after the target site becomes empty. Let us see the example in Fig. 3. 
The configuration does not change until the first particle moves to the site one. It took 1 
step before the first move and hence wi^i = 1. Now two particles can move to the right. 
It took 1 step (resp. 2 steps) before the second (resp. first) particle makes the first (resp. 
second) jump and hence wi^2 = 1 (resp. W2,i = 2). After the second jump of the first 
particle the second particle waits 2 steps before the next move and hence ^2,2 = 2. One 
continues this process to get all Wjj's for a given sample of the dynamics of dTASEP. For 
this example the matrix looks like 



w 



1113 
2 2 10 
10 1 
12 10 



(3.1) 



It is clear that the time evolution of the dTASEP is equivalent to the infinite matrix Wij. 
If we call a matrix whose elements are nonnegative integers an N matrix, is an N matrix. 
When considering various realizations of dynamics, each Wij is a random variable and 
geometrically distributed with parameter g; F{wij = /) = (1 — g)g' (/ = 0, 1, 2, ...). 
Next let us define G(M, A^) by 



GiM,N)= max {J2 ^ij}- (3-2) 



Here vr is a path from (1,1) to (M, A^) which consists of elementary path connecting points 
+ or — > + 1). In the matrix representation as in (13.11) . vr is a 
down/right path, i.e., it starts from (1,1) and goes either down or right until it reaches 
{M,N). IIm,n is the collection of all such paths. If we regard Wij as representing a 
value of a potential energy at a position (13. 2p can be considered as the problem of 

the zero temperature directed polymer in random medium. One can confirm oneself that 
G{M, N) + M + A^ — 1 is equal to the time at which the A^th particle has just made the 
Mth step. Considering the dynamics of the A^th particle, one notices that A^(t) > A^ is 
equivalent to saying that the A^th particle makes the A^th jump before t. Hence we have 

F[N{t) >N]= ¥[G{N, N) + 2N-l<t]. (3.3) 
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Though w is in principle an infinite matrix, to study N{t), we can restrict our attention 
to the N X N submatrix {wij}i<ij<N of w. This is related to the Markov property and 
the fact that in the TASEP each particle can not affect the dynamics of the particles on 
its right. For instance, from the 3x3 submatrix of w, 



1 1 2 

2 2 1 
1 



(3.4) 



one sees that G(3, 3) = 6 with the maximizing path, (1, 1) — » (2, 1) — > (2, 2) (2, 3) — > 
(3, 3). This gives G(3, 3) + 5 = 11 which agrees with the time at which the third particle 
has made the third hop in Fig. 3. 

To proceed we use a bijection known as the Robinson-Schensted-Knuth(RSK) algorithm 
between an N matrix and a pair {P,Q) of semistandard Young tableaux(SSYT). A short 
explanation about the SSYT and the Schur function is given in the appendix. For more 
details and the RSK algorithm, see for instance [19,20]. The pair of SSYTs corresponding 
to (I33D is 
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In the language of SSYT, the G{N, N) is given by 

G{N,N) = Ai 



(3.5) 



where Ai is the length of the first row of P. This is easily seen as a property of the RSK 
algorithm. For our example, Ai = 6 is the same as ^(3, 3). Therefore the problem of N{t) 
is reduced to a combinatorial problem of SSYTs. If we denote the number of SSYT of 
shape A and entries from {1, 2, ... , A^} by L{X, N), we have 



F[G{N, N) + 2N-l<t] 



X:\i<t-2N+1 



gi^isA(l,...,l,0,...) 




(3.6) 



where a formula flA.4p is used in the second equality. Rewriting (13.61) using the formula 
flA.6p and hj = Xi — j + 1, one obtains 



^ t-2N+l 

F[G{N, N) + 2N-l<t] = - 



Z ^ U 

hj=-N+l l<j<l<N 



N 



(3.7) 



Here and in the following the symbol Z is used to represent a normalization constant. For 
the present case it is 

i7V(7V-l) N 

Z = ^^ ^U^Kj-^V- (3.8) 
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The expression of the RHS of (13 .yp has a strong similarity with a quantity appearing in the 
random matrix theory. Hence generahzing and applying techniques from random matrix 
theory, one can study the asymptotic behavior of the quantity to obtain 



lim P 



N- 



N 



< s 



(3.9) 



with d = g^/^(l + y/qy^^{l — ^) ^. In terms of N(t), this is equivalent to 



hm P 



N{t) 



> 



(3.10) 



with d' = 2-^/3gi/^(l - g)^/^ r^j^g function Fa appearing on the RHS is called the GUE 
Tracy- Widom distribution function [21] and describes the distribution of the largest eigen- 
value of the Gaussian unitary ensemble (GUE). Hence the above equation states that the 
current of the TASEP for the step initial condition is, in the appropriate scaling limit, 
equivalent to that of the largest eigenvalue of the GUE. It is stressed that F2{s) is com- 
pletely different from the error function which is the distribution function for the ordinary 
Gaussian fluctuation. This is the most basic and fundamental result in the approach of 
random matrix theory to the ASEP. Basics of random matrix techniques to obtain F2 and 
related quantities will be given in the next section. 

Before ending the section, let us consider the cTASEP limit. In this limit, Wij becomes 
an exponential random variable with parameter one and (13. 7p reduces to 



F[GiN,N)<t] = - 



N 



n {xj-xiyY[e-''^d''x 

l<j<l<N j=l 



(3.11) 



with 



N 



Z = l[j\{j -1)1 (3.12) 

i=i 

The same quantity appears in the random matrix theory as the distribution of the largest 
eigenvalue in an ensemble called the Laguerre unitary ensemble. As a result, for the 
cTASEP, we have 

-N{t)- 



lim P 

t— »oo 



2-4/3^1/3 



> -s 



F2{s). 



(3.13) 



4 Random matrix techniques 

In this section, we explain the basics of the techniques from random matrix theory which 
would be helpful to understand the results for the ASEP. In physics the random matrices 
were first introduced to model a complex Hamiltonian of nuclei. But the random matrix 
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theory has also found a lot of applications not only in nuclear physics but also in other 
branches of physics, mathematics, statistics and so on. The basic reference of random 
matrix theory is [22]. The book by Forrester [23] is also very useful. 

A random matrix is a matrix with random elements. There are infinitely many variety 
of random matrix ensembles but there are some random matrix ensembles with especially 
good properties. Among them are the Gaussian ensembles which are defined by the prob- 
ability distribution of N x N matrices, 

P{H)dH = e-^^^'dH. (4.1) 

Here f3 takes three different values (3 = 1,2,4. When (3 = 1, H is taken to be a real 
symmetric matrix and dH = Y[f=i ^^jj Y[j<i dHji. Writing down the Tr in (14. ip explicitly, 
one sees that this definition is equivalent to saying that each independent element obeys a 
Gaussian distribution with variance 1 for the diagonal and 1/2 for the off-diagonal elements. 
This ensemble is called the Gaussian orthogonal ensemble(GOE). For j3 = 2, H is taken to 
be a Hermitian matrix and dH = Y[f=i ^^jj Y[j<i ^^ji^^ji 'w^here H^^ and Hji are the real 
and imaginary parts of Hji respectively. This is the Gaussian unitary ensemble (GUE). 
For /5 = 4, the ensemble is called the Gaussian symplectic ensemble (GSE). For this case, 
H is taken to be a self-dual quaternion matrix, whose definition is not given here. 

In many applications we are interested in the statistics of eigenvalues. For the Gaussian 
ensembles, we can derive the joint density function for eigenvalues explicitly. Basically 
we change the coordinate from the original matrix elements to the eigenvalues and the 
diagonalizing matrix and then integrate out the latter degrees of freedom. The procedure 
is explained in [22,23] but is omitted here since it is not directly used in the applications 
to ASEP. The resulting joint probability density function of eigenvalues {xj}j=i^...^iv is 

1 ^ 
PM^,,...,x^) = - H Ixj-xil^He-f'^V^ (4.2) 

l<j<l<N j=l 

for (3 = 1,2,4. Let Fj^f^ denote the corresponding measure. Using this expression one can 
study in principle any statistical property of eigenvalues of the Gaussian ensembles. Here 
we focus on the statistical behavior of the largest eigenvalue. From the joint eigenvalue 
distribution (14.21) . the probability that the largest eigenvalue xi is less than u is given by 

1 r ^ 
^N^x, <u] = - U - "^il^ n e-^'^'^^'dx,. (4.3) 

Next for the case of /5 = 2 we rewrite this expression to another which is more suitable 
for asymptotic analysis. The probability density (14.21) for (3 = 2 can be written as 

i det(y.,(xO)S=, det(^,(xO)S=i (4.4) 

if we notice 

n {xi-x,)=Aei{x^-Xi=i. (4.5) 

l<j<l<N 
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and set v^j(a;) = i'jix) = ^'^ 1'^ . We show the identity 



- / c?xdet(v.,(xO)S=idet(^,(xz))S=in(l + ^(^^)) = ^^*(l + ^^)- (4.6) 
Here the determinant on the RHS is the Fredholm determinant defined as 

°° / -\\m poo poo 

det(l + Kg) = Y^ ^ / dxi--- dx^g{xi) ■ ■ ■ g{xm) det {K{xj,xi))]^i^^, (4-7) 

where the kernel is 

AT 

K{xi,X2) = i'jixiM'\mix2) (4.8) 

j,l=l 

with the matrix A being 

POO 

Aj^i = / dx(pj{x)ilJi{x). (4.9) 



There are several ways of proof; here we follow [24]. Using the Heine identity, 

dxi ■ ■ ■ J dxN det{(fj{xi))fi^idet{ipj{xi))fi^i = det( J dx(^j{x)ipi{x))^i^i, (4.10) 
the LHS of (14.61) can be rewritten as 

det{Aj i + J dxLpj{x)'ipi{x)g{x))/ det A 
= det{lji + j dx{A-^^)j{x)g{x) ■ ^/(x)), (4.11) 

where 1 is the identity matrix of size A^. One can regard this as det(l + BC) if we set 

Bj,, = {A-^^),{x)g{x), a,, = ^j{x). (4.12) 

Then using a simple identity det(l + i?C) = det(l + Ci?) and rewriting B, C again in terms 
of (p,ip, we get (14.61) . Some manipulations here are rather formal but can be justified [24]. 
Since (14. 3 p with /3 = 2 is written in the form of LHS of (14.61) with g{x) = —X{u,oo){x) where 
X{u,oo){x) = l(x > u),0{x < u), the probability that the largest eigenvalue Xi is less than 
u in the GUE is written in the form of the Fredholm determinant, 

FN2[xi<u] = det{l + Kg). (4.13) 

To proceed, a nontrivial task is the computation of the inverse matrix A~^. For the 
GUE, this is (implicitly) done as follows. Let us remember that the addition of a row or a 
column to another does not change the value of the determinant. Then the functions ipj 
and Ipj in (14.41) . which we have taken to be the monic polynomials x^~^, can be replaced 
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by any polynomials of the same degree modulo a numerical factor. For the GUE, it is 
convenient to take 

^,{x) = ^|J,{x) = e-^'/'H,{x) (4.14) 
where Hn{x) is the Hermite polynomials [25], 

H^{x) = {-ire^'^e-^\ (4.15) 
Using the orthogonality relation of the Hermite polynomials, one sees that A is diagonal, 

/oo 
Hj{x)Hi{x)e-''^dx = VT^2^j\5ji. (4.16) 
-oo 

Then the kernel (14. 8 p for the GUE, which we denote by Kn2{xi,X2), is calculated to be 



N-l 



?i=0 



This expression of the kernel allows one to study various statistical properties of the eigen- 
values such as the eigenvalue density, the two point correlation function, etc. To take a 
scaling limit for the largest eigenvalue, we use the asymptotic expression of the Hermite 
polynomial [26], 



e-5-'i/„(x) ~ 7rh(2"+i)/^(n!)i/2n-i/i2Ai(2i/2^i/6(x - v^)), (4.18 
where Ai(x) is the Airy function defined by 

nco+it 

27r 



Ai(a;) = — / d'^'+i'^dz (4.19) 



— oo+ie 



with e > 0. If we set Xj = V2N + ^j/{V2N'^/^) (j = 1, 2), one has for large 

Kn2{xi, X2) ~ y2Ari/6/C2(ei, 6) (4.20) 

where 

r-(,,^ Ai(ei)Ai-(6) - Ai(6)Ai-(6) 

/C2(4i,6) = 7 7 (4.21) 

is known as the Airy kernel. Hence one obtains 



lim 



(xi - V2N)V2N^/^ <s\= det(l + /Ca^) =: ^2(5). (4.22) 



The function F2{s) is defined as the Fredholm determinant (14.71) with the kernel /C2 (I4.19P 
and G{x) = — X(s,oo). The function -^2(3) describes the fluctuation of the largest eigenvalue 
of GUE in the scaling limit and is called the GUE Tracy- Widom distribution function. It is 
known to have a different representation using the solution to the Painleve II equation [21]. 
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Equation (14.221) says that the largest eigenvalue xi in the GUE behaves like xi ~ V2N in 
average, has a fluctuation of order 0{N~^^^) and has the F2 as the distribution function 
in the N 00 limit if appropriately scaled. 

The function F2{s) is the same as the one which appeared in the RHS of (13. 9p and (13.101) . 
Notice a close similarity between (13. 7p and (14.31) with j3 = 2. The main difference in (13. 7p 
is that the Gaussian {x G M) is replaced by {h = 0, 1,2, . . .). As a consequence 
we use the orthogonal polynomials which are orthogonal with respect to the weight q^, 
which is known as the Meixner polynomials [27]. Then (13. 7p could be called the Meixner 
ensemble representation for the current of the TASEP. The derivation of (13. 9p from (13. 7p 
is more involved but is similar to the one for the GUE [18]. 

For GOE(/? = 1) and GSE(/3 = 4) the power of the determinant in (14. 3 p is different 
from 2 and the above calculation does not work. But it is possible to modify the arguments 
and obtain analogous results [24, 28] . They read 

= Ff^is) (4.23) 

with /3 = 1,4. The functions Fi{s) (resp. -^4(5)) is called the GOE (resp. GSE) Tracy- 
Widom distribution function. They are expressed as Fredholm determinants with 2x2 
matrix kernel. The Painleve II representation is also known [28]. 

A remark is in order. To compare the results of ASEP with those of the Gaussian 
ensembles, it is sometimes more convenient to change the normalization from (14. ip to 

P{H)dH = e-^^'^'dH. (4.24) 
Then Km and Km are replaced by 

= Fpis) (4.25) 

for /3 = 1, 2, 4 in which the exponents of are the same as in (13.90 . 

Let us now consider a time dependent version of the random matrix ensemble known 
as the Dyson's Brownian motion [29]. The random matrix H{t) now depends on time t. 
Each independent component of the matrix performs the Ornstein-Uhlenbeck process for 
which the Fokker-Plank equation for the probability density p{t, x) takes the form, 

^^^(^' ^) = ^) + ^))- (4.26) 

Here D is the diffusion constant and should be taken appropriately depending on the 
ensembles in question. Let us consider the stochastic dynamics of the eigenvalues of this 
random matrix. We are interested in the multi-time joint distribution of the eigenvalues. 
It is known that, for general (3, time evolution of the probability density of the eigenvalues 
is described by an equation which can be mapped to the quantum Calogero model [23]. 
For the GUE(/9 = 2) case, the system is a free fermion and we can get detailed information 



lim 



P 



Af/3 



[Xi 



lim P 



N(3 



xi - 2N 
A^i/3 



< s 
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from this fact as we see below. For the G0E(/3 = 1) and the GSE(/? = 4), however, the 
analysis becomes much more involved and the computation of correlation functions have 
so far been very restricted. 

Let us consider the stationary situation for the GUE in which the dynamics and the 
initial conditions are both taken to be the GUE type. We call this the tGUE (This 
is not a standard notation). If we denote the eigenvalues of the matrix Hji{t) by x* 
(j = 1, 2, ■ ■ ■ , A^), the dynamics of them look like Fig. 4. The joint density function at two 
times can be written product of determinants, 

(4.27) 

where ^f'\x) = ^f'\x) = e-^'/^^^-i = 2, . . . , N) and 



(pt„t2{xi,X2) = -y======exp 

a/(1 - e-2(t2-*i))7r 



{X2 - e-(*2-ti)xi)2 
1 _ e-2{t2-ti) 



(4.28) 



This is a generalization of (14.41) to the two time case. The derivation can be found in 
[22, 23] but is omitted here since it is again not directly used in the applications to ASEP. 
Furthermore the joint probability that the largest eigenvalue Xi (the height of the top 
curve at tj in Fig. 4) is less than uj at time tj [j = 1,2) is again given by a Fredholm 
determinant. Let us introduce a matrix 

/OO /"OO 
dxi / dx2(fJ^\xi)(f)ti,t2{xi,X2)i'f''\x2), (4.29) 
■OO J ~oo 

and 

/OO 
dxi(ff^\xi)(f)t^^t2{xi,X2), (4.30) 
-OO 
POO 

^Pp\x,)= / dX2(l)t,Mi^l,X2)^f'\x2). (4.31) 



Then we have 



^N2[x{' < ui, xf < U2] = det(l + KN2g) (4.32) 



where the kernel is 

KN2{ti,xi;t2,X2) = KN2{ti,xi;t2,X2) - 0ti,t2(a;i,a;2), (4.33) 



N 



KN2{ti,xut2,X2) = J2^\''\^iM~']jM'"\x2), (4.34) 

g{tj;x) = -X(u,,oo){x), (j = l,2), (4.35) 

with the convention that (j)t-^^t2{xi, X2) = when ti > t2. In this kernel (^1,^2) should 
be understood as a representative of either (ti,ti), (^1,^2), (^2,^1), (^^2,^2)- Similar abuse of 
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notation will be used in rest of the paper as well. The formula (14.321) can be derived by a 
generalization of the argument to get (14. 6p . 

For the tGUE, the matrix A in (I4.29p is again diagonalized by taking (pj^^ and ^/'j*^'' to 
be the Hermite polynomials. Then the kernel turns out to be 



KMtu xi; h, X,) = e-^(^?+^i) ^^^^^i^^e-("+V2)(..-.0 (4.36) 



n=0 



Using the integral formulas for the Hermite polynomial [25], we can rewrite this in a double 
integral form, 



KN2{ti,Xi;t2,X2) = e y) _— _ / dzi / dz2 — 

[Zin) J Jy Z2 



N 2-^6*1 -*2 - Z2 

(4.37) 

where F is a contour enclosing the origin anticlockwise and 7 is an any path from —ioo 
to ioo. This representation of the kernel has a form suitable for the saddle point method. 
Performing the asymptotic analysis, one gets with tj = 



lim Fm 



(x^^ - Vm^ V2iV^/^ < Sj {] = 1,2) = det(l + /Cs^). (4.38) 
Here Q{Tj;x) = — X(sj,oo)(a;) (j = 1,2) and the kernel is 

/C2(ri, Ci; r2, 6) = ^2(ri, 6; r2, 6) " Mn, 6; rs, 6) (4.39) 

where 

<;Ae-^l"-''>Ai(6 + A)AiK2 + A), (4.40) 

*.(..«.;..6)4P'''"'"""'*'«'"'>*"«^^'>- ^T^- (4.41) 

[0, Ti > T2. 

This is a generalization of (I4.22p . This describes the stochastic dynamics of the largest 
eigenvalue of the tGUE in the scaling limit and is called the Airy process [12,30]. 

The weight (14.270 is for the stationary situation. One can also consider the case where 
the dynamics is still (3 = 2 but the initial condition is taken to be GOE or GSE. This 
describes the transition from the GOE or GSE to the GUE. For these transition ensembles, 
a similar analysis as above is possible [31]. For the dynamics of the G0E(/5 = 1) and 
GSE(/5 = 4) cases, however, much less is known as already mentioned. 



5 Discrete polynuclear growth model 

In section 3, we have obtained a formula for the current fluctuation by mapping the problem 
to that of combinatorics of Young tableaux. In this section we derive a generalized result 
by introducing the discrete polynuclear growth (dPNG) model [30]. 
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The dPNG model can be defined by interpreting G{i,j) as a height of a surface. Let 
us introduce a new coordinate (Fig. 5), 

s = i+j — 1, x = i— j (5.1) 

and set h^{s) = G{i,j). In the dPNG model s and x play the role of time and space 
respectively and h^{s) is considered as the height of a surface at time s and at position x. 
From the definition of G{M, N) (13.21) . the growth dynamics of h^{s) is given by 

h%s + 1) = max{/i^~^(s), h'^+^s)} + (5.2) 

Here Wij is the waiting time matrix w = {wij} and is taken to be zero if (s + x)/2 + 1 or 
{s — x) /2 + 1 is not an integer. The time evolution of h^{s) until just before s = 3 for our 
example is depicted in Fig. 6. 

The above definition can be restated as follows. At time s = 0, the surface is flat; 
/i^'(O) = for all X. At each time s, nucleations occur at sites x = —s + 1, —s + 3, . . . , s — 
3, s — 1. A nucleation at each site is independent and the height of a nucleation at each 
site obeys the geometric distribution with parameter q. In the meantime the surface also 
grow laterally in both directions deterministically with unit speed. Sometimes two steps 
collide whence the height is taken to be the higher one. The dPNG model is a model of 
stochastically growing surface. 

We further introduce the multi-layer version of the dPNG model as follows. Suppose 
that there are infinitely many layers below the original surface. We denote the height of 
the jth layer by hj{s) {j = 1, 2, . . .). In particular h^{s) = h^{s). Initially they are equally 
spaced as hj{0) = 1 — j. The growth of the j'th [j > 2) layer is determined by the layer 
above. When there occurs a collision with height I at {j — l)th layer, then there occurs a 
nucleation with height / at the same position of the jth layer. They also grows laterally 
in both directions with unit speed but there are no stochastic nucleations. This is the 
multi-layer PNG model. A snapshot of a Monte-Carlo simulation is shown in Fig. 7. 

For a given time s if we regard the PNG space coordinate x as time and the height 
coordinate h as space, the multilayers can be considered as a system of random walkers 
with the non-colliding condition. This kind of object has been known as vicious walks 
in statistical physics [32]. The weight of the non-colliding walkers (or, line ensembles) 
is known to be simply given by a product of determinants with entries being transition 
weight of each walker. Mathematically this follows from the Lindstrom-Gessel-Viennot 
theorem [33,34]. Since our main interest is to study N{t) and h^{2N — 1) = G{N,N), let 
us fix s = 2N — 1 and write = h^{2N — 1) in the sequel. The transition weight of each 
walker at each time step is given by 

<Px,x+i{hi,h2) = < ^ ^, (5.3) 



when X is odd and 



h2 < hi, 
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when X is even. This is a consequence of the fact that a nucleation of a height / occurs 
with probabihty proportional to qK The transition weight from Xi to X2 is given by the 
convolution from through (f)x2-i,x2 is compactly given by an integral formula, 



0, Xi> X2 



1 > 2 
(5.5) 

where Cr is the contour enclosing the origin anticlockwise with radius R. The weight of 
the vicious walks that their positions are {hj} is simply given by a product of determinants 
of this transition weight, 

2N~2 

n det(0f,+,(/.J,/.r')) (5.6) 

x=-2N+l 

with the initial and final condition hj'^'^^^ = h'^'^~^ = 1— j {j = 1,2,...). Since the 
number of layers is infinite, the determinant on RHS is in principle infinite dimensional. 
But it can be regularized by restricting the number of layers to be finite for a moment and 
then taking the limit. 

Notice a similarity of this weight (15. 6p to the joint density of tGUE (14.271) . One can 
use a similar argument to express the joint distribution of the top layer in the form of the 
Fredholm determinant. It reads 

P[/if < ui, /if < U2] = det(l + K^^g) (5.7) 

where the kernel is (for even xi,X2) 

ir^P(xi, X2, h2) = kf{x^, h; X2, /^s) - 0^L,(/ii, h2) (5.8) 



with 



Here 



K^''{x^,h;x2,h2) = J2^j"\h,)[A-'],M^'\h2). (5.9) 



^.V = &+i,2^-i(l-J>l-0, (5.10) 

^t\(/^) = 0lWM(l-J», (5-11) 
^f-\W = C^-i(/^,l-j), (5.12) 

for j,l = 1,2, . . .. The next step is the computation of the inverse matrix A~^. For the 
present case, it can be obtained by utilizing the fact that A is a Toeplitz matrix, leading 
to the kernel [18], 



Kff{xi,hi;x2,h2) = / ^/ dz; 



7g)^2-^i f dzi f 4^-^ (1 - yg/zi)^+^i/2^1 - ygza)^"^^/^ 



{2my Jc,^ 4^ Jc,^ ^''z,- Z2 (1 - Vg;.i)^--/2(1 - yg/z^)^^-/^ ' 

(5.13) 
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where ^ < R2 < Ri < l/^ [30]. This is similar to the kernel for the tGUE 

Our original interest was to consider h^{2N — 1) which is directly related to N{t). But 
in the surface growth picture it is natural to extend our analysis to the equal-time (in s) 
multi-point distribution because it is related to the roughness of the surface. One can show 
that the scaled height 

h^{s = 2N-l)-a{bo+'-^)N 



with X = 2boN + 2c(6o)-^^^'^t tends to the Airy process as N 00 where 



(5.14) 



a{b) = ^ + v^r^^) , (5.15) 



2^ 

Q 
1 

d(b) = — -(Vl + b + Ml - 6)5(Vl - & + Ml + &)^ (5.16) 

(l-g)(l-62)r 

c(6) =q-'^{\-b^f^{^fYTb^ ^^^Y^)^{^fY^^ ^^fYTb)h. (5.17) 

For 60 = this was proved in [30]; the general b case is obtained as a special case of the 
results in [35]. Also notice that if we take x = and consider the fluctuation at the origin, 
this reduces to (13. 9p . It is possible to translate this result into the language of the TASEP 
by noticing that /i^(s) corresponds to the time at which the ^"^"'"^ th particle has made 
the ^i|^th hop. It gives a joint distribution of the times at which ASEP particles pass 
certain bonds. It should be remarked here that the equal-time in PNG picture does not 
correspond to the equal-time in TASEP. 

In our argument the multi-layers were introduced as an auxiliary object which makes 
it possible to analyze the properties of the original surface related to the quantities in 
TASEP. From the random matrix point of view, the multi-layers have a meaning that in 
the limit the behaviors of the jth layer is the same as the jth largest eigenvalue of GUE. 
Moreover if we set 

h] = \]-3 + \ (j = l,2,...), (5.18) 

the probability distribution of {A^} is known to be a special case of the Schur process [36,37] 
which assigns the weight of the form, 

Sa(i) (Po )Sa(i) //.(I) (Pr)5A(2)/M(i) {P\) ■ ■ ■ S\W/f,(N^i) ipN_i)Sx(N) (p^), (5.19) 

on the set of Young diagrams satisfying 

A(^) D /i(^) C A(2) D ■ ■ ■ D C A(^). (5.20) 

Here Sp/^(p) is the Schur function with the specialization of algebra p. In particular for 
X = 0, the probability distribution of A = {A°}j=i 2,... is proportional to 

SA(Vg,...,v^,0,...)^ (5.21) 

N 
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Restriction of our attention only to Ai gives another derivation of (13. 6p . 

The hmiting case of cTASEP has been studied in detail in [38,39]. 

In the PNG picture, there is another limit of much interest: q ^ 0, N —>■ oo with 
■y/qN fixed. This corresponds to the continuous PNG model in which time s and space x 
coordinate are continuous and the height h takes integer values. In the study of the ID KPZ 
systems using random matrix techniques, this model has been studied most extensively 
[11,12,40-42]. From the universality there is little doubt that the continuous PNG model 
and the TASEP give the same results for the quantities which are not model dependent 
and characteristic of the KPZ universality class. 

6 Variants 

So far we have concentrated on the current fluctuation for the step initial condition. It is 
also possible to generalize the analysis to some other initial conditions, boundary conditions 
and quantities. 

6.1 External sources 

One can generalize initial conditions to the P-,p+ situation in which the sites on the 
negative side are occupied with the stationary measure with density p_ and the sites on 
the positive side with the one with density p+. Of course p_ = l,p+ = corresponds 
to the step initial condition. In the stationary measure with density p the waiting time 
of each particle is geometrically distributed with parameter 1 — J(p)/p. In the language 
of G{i,j) this situation can be realized by setting G(l,j) (resp. G{i,l)) as a geometric 
random variable with parameter 1 — J{p^)/ p_ (resp. 1 — J{p+) / p+). This interpretation is 
explained in more detail for the cTASEP case in [41]. In the dPNG model this corresponds 
to the model with external sources. The fluctuation of the height at the origin was studied 
in [43] and those at multi-points in [44]. They are readily translated into the language of 
the ASEP. 

6.2 Stationary Two point function 

One of the most important achievements in the application of the methods of random 
matrix theory is that the stationary two point function for the ID KPZ universality class 
has been computed quite explicitly by Ferrari, Prahofer and Spohn. In the PNG picture 
this is related to the limiting case of p_ = p+ for the model with external sources above. 
We do not go into details here. See the references [39,41,42]. 

6.3 Alternating initial condition. 

Another tractable initial condition is the alternating initial condition in which initially all 
even sites are occupied and all odd sites are empty (Fig. 8). In the dPNG picture the surface 



16 



is initially flat. The nucleations occur at all even sites at even time and at all odd sites at 
odd time. This should be contrasted to the step case in which the nucleations at time s 
occur only at sites \x\ < s. In addition although the alternating case macroscopically looks 
the same as for the stationary case with density p = 1/2, it turns out that the fluctuation 
is different. 

The fluctuation of G{N, N) can be studied by using a symmetry property of the com- 
binatorial approach and is shown to be Fi{s) which appeared in fl4.25p as the largest 
eigenvalue distribution of the GOE [45,46]. It is also possible to define the multi-layer 
PNG model. From the symmetry, the weight of A defined from the height of multilayers 
by f lS.lSp is given by 

sa(v^,...,v^,0,...) (6.1) 

N 

where partition A is restricted to ones such that it has the form A = 2/x where n is another 
partition. From this one can show that the fluctuations of the multi-layers at the origin 
are the same as the largest eigenvalues of GOE. Combined with the fact that the flat case 
is statistically translationaly invariant, this suggests that the multi-layers are the same as 
the GOE Dyson's Brownian motion (tGOE) as conjectured in [47]. But unfortunately the 
analysis of the multi-layers for the flat case has been difficult. For the multi point joint 
distribution, at least for now, we need to employ a little different method which will be 
explained in subsection 7.4. 

6.4 Half-infinite system 

In the study of the ASEP, the boundary effects are important. The random matrix tech- 
niques also allow to study some properties of the half-infinite system in which the particle 
can enter the system at the origin with probability ■j^yq (0 < 7 < l/^/q) if the site is 
empty. In the language of G{i,j) this can be realized by putting the symmetry condi- 
tion G{i,j) = G{j,i). The symmetrized version of combinatorial approach was studied 
in [45,46]. As the fluctuation of G{N, N) there appear Fi{s) and F^i^s) in fl4.25p depending 
on the value of 7. In the ASEP language this says that the fluctuation of the current at 
the origin, i.e., the number of particles which entered the system in time t is equivalent to 
that of the largest eigenvalue of the GOE and GSE in the scaling limit. The current at 
other points were studied in [48]. There the fluctuation of the current near the origin is 
shown to be described by the transition ensembles mentioned in section 3. 

6.5 Position of particles 

So far we have mainly studied the fluctuation of the current, i.e. the number of particles 
which crossed a given bond until time t. It is also interesting to consider a related quantity, 
the position of particles at a given time t. In this section, we consider the step initial 
condition. 

Let us denote the position of the iVth particle at time t by XAr(t). If we remember that 
G{M, N) + M + — 1 is the time at which the A^th particle has made the Mth jump. 



17 



properties of XN(t) for a single particle can be obtained by a translation of results for the 
current. If we take t ^ oo limit for a fixed N, it is known that XN{t) goes to the dynamics 
of the largest eigenvalue of GUE [49]. Here we consider another limit t oo, N oo with 
u = N/t fixed. The quantity x^t{t)/t has a deterministic limit, 



lim{x^t{t) +ut-l)/t = a*{u) (6.2) 

t— >oo 



where 



^,^^^^\l-q-il-2q)u -2^1^(1- iy)q{l-q), < < 1, 
1 0, u > 1. 

For a single particle, the fluctuation of XAr(t) around this position can also be obtained 
from the result for the current by setting r = in fl5.14p . 

It is further possible to study the equal time joint distribution of the position of several 
particles by introducing a multi-layer PNG model which is a little different from the one in 
section 5. See Fig. 9. It is defined as follows. Let k^(t) be a height of the jth (j = 1, 2, . . .) 
layer at time t (= 0, 1,2,.. .) and at position v = 0, 1, 2, . . .. Initially they are all flat; 
kj{0) = 1 — j for all v. At each time t at each layer nucleations of height one may occur at 
sites f = l,3,...,2t — 1 with the condition that the layers do not touch. Here a nucleation 
at each site is independent, has a shape depicted as in Fig. 9, and occurs with probability 
q. In the meantime it grows laterally to the right with speed 2. Notice a similarity and 
difference between the PNG model here and that in section 5. This is the same as the 
DR-paths in the Aztec diamond [50,51]. One can confirm oneself that the dynamics of the 
top layer k^(t) = kKt) gives, if set 

XN{t) =t-2N + 2-k'^^~\t), (6.4) 

the correct statistical properties of the TASEP particle positions. 

The multi-layers in this PNG model can again be regarded as vicious walks. The 
transition weight of each walker is now 

<Pv,v+iiki,k2) = < (6.5) 
[0, fca < fci, 

with p = q/{l — q) when v is even and 

€,.+iih,k2) = l' ' (6.6) 

10, otherwise, 

when V is odd. The weight of the vicious walks are simply given by a product of determinant 
of this transition weight, 

2Ar-l 

n det(</);,+,(fcJ,fcr^)) (6.7) 

v=Q 
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with the initial and final condition = k'j^ = 1 — j (j = 1, 2, . . .). 
The joint distribution has the form, 

¥[xr,it) >t-2ri + 2-ui,Xr2{t) >t-2r2 + 2-U2\ 
= F[k^''-\t) < ui, k^'-'-^t) < U2] = det(l + K^g) (6.8) 

where the kernel is given by 

K*{ri,ki;r2,k2) = ^^(ri, /ci; rs, /C2) - ^^i.ral^i' ^2), (6.9) 

with 

iff,, ^2 {i + i/ziy~'^+\i-pz,Y^ 



Iff z 

K*{ri,ki;r2,k2) = ——T / dzidz2—^ 



4^+' ^2 -^1(1 + l/z2y-^'^+\l - VZ2Y- ' 

(6.10) 



cj^.^2)-^y^^c^.^^^^) (6.11) 



where Ri < R2 and < i? < 1. Applying the asymptotic analysis, one can show that 

a*{h)t-ht + l-XN{t) 



with N = ut + 2c*(z/)t^/^r tends to the Airy process where 



(6.12) 



a'{i^) = (1 - - (1 - 2qy - 2-^v(\ - v]q(\ - q], (6,13) 



<,» = (.(.-i))M(i-#|n-,/^ V^-y^l . (6") 



.(.,.(.(.- D). I l.^^) (;ii^-y^ I . (6,15, 

In the cTASEP limit, (16321) with (l635ll replaced by 

a*(z/) = (v^-l)^ (6.16) 
d\v) = u'/%^/T/^ - 1)2/3, (6.17) 
c*{u) = ul{^/l/^-iy/^ (6.18) 

tends to the Airy process. 

What have been explained in this subsection is the equal-time multi-particle joint dis- 
tribution. It is also possible to study the multi-time joint distribution of the position of a 
single particle (tagged particle) [49]. 



19 



7 Green's function method 



So far we have explained how the fluctuation is computed by combining a combinatorial 
argument and techniques from random matrix theory. There is a related but a somehow 
different method using the Green's function for the TASEP. This is the probability that 
the particles starting from sites ...,?/ at {yi > 1/2 > ... > Hn) at time are on sites 
xi, . . . ,xn (yXi > X2 > ... > xn) at time t and is also called the transition probability. We 
denote this quantity by ^(xi, . . . , xn] t\yi, . . . , z/n; 0) or, when we fix the initial condition, 
by G{xi, . . . ,XN',t). For the cTASEP, a determinantal formula was obtained by Schiitz [52]. 
In this section we explain some applications of the formula. In this section we only consider 
the continuous time case. For the discrete case the expression for the Green's function was 
obtained in [53] but has not been exploited to obtain the results in the scaling limit. See 
also [54]. 

7. 1 Formula 

First let us remember the formula [52]. It reads 

G{xi,X2,--- ,XN;t) = det[Fi_j{xN-i+i - yN-j+i;t)]j,i=i,-,N (7.1) 

where the function Fn{x,t) appearing as matrix elements of the determinant is explicitly 
given by 

F.X.-.t)^e-'L^yf^L. (7.2) 

To see the structure we write down the determinants for = 1,2. For A^ = 1 this is 
simply 

G{xi; 0) = Foix ~y;t)= / ' e"*. (7.3) 

{x - yy. 



For N = 2 the formula gives 



G{xi,X2;t) 



Fo{x2-y2;t) Fi{xi-y2]t) 
F_i(x2-yi;t) Fo(xi-?/i;t) 



(7.4) 



The proof of (17. ip consists of checking the master equation and the initial condition which 
should be satisfied by the Green's function. 

Some formulas of the function Fn{x]t) which are useful in the following are given. It 
has an integral representation, 

FJx;t) = / dw- r^-re"* (7.5) 

' 27iiJ^^_^ (l + «;)"+-+i ^ ^ 
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where Fq,-! represents a contour enclosing both w = and w = —1 anticlockwise. They 
satisfy the recurrence relations, 

j^Fnix; t) = Fn{x - 1; t) - t), (7.6) 

/ rftF„_i(x-l;t) = F„(x;t), (7.7) 
Jo 

oo 

F4a;;t) = ^F„_i(a;i;t), (7.8) 

X\=X 

F„_i(x; t) = F„(x; t) - F„(x + 1; t). (7.9) 
7.2 Position fluctuation of a single particle 

Let us first study the fluctuation of the position of a particle. A basic observation is that 
the probability we are interested in is written as a sum of the Green's function, 

f>[xN{t)-yN>M]= J2 G{xu--- ,XN;t\yi,--- ,yN;0). (7.10) 

j/]V+Af<a::]v<---<a;i 

Here P indicates the measure of the TASEP for our current initial condition and the 
summation on RHS is over all x^'s which satisfy the constraint Vn + M < xn < ■■■ < xi. 
The next thing to do is to rewrite the RHS to a form which is suitable for an asymptotic 
analysis. In [55,56], it was shown that the quantity can be written in the form of a multiple 
integral which reduces to (13.111) for the step initial condition. 

In this subsection we state a similar result for the fluctuation of x^it) for the step 
initial condition, i.e., ior yj = 1 — j(j = 1, . . . , A^). We have 

oo oo N X. 

F[x^{t)+N-i>M]=-j2---T. n (^^■-^')'n^ (7-11) 

xi=M xj^=M l<j<l<N j=l ^' 

with 

TV 

^^^iV(iv-i)/2-Q^.,_ ^7^^2) 

Since the orthogonal polynomials for the weight * is called the Charlier orthogonal 
polynomials [27], this expression could be called the Charlier ensemble representation of 
the fluctuation of X]\f{t). The proof is almost the same as in [55] and is omitted here. But 
to illustrate the method we give a step-by-step calculation for the N = 2 case. It proceeds 
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as follows. 



P[x2(t) > M- 1] = Yl G{xi,X2;t\0,-l-0) 

AI-l<X2<xi 

Fo{x2 + l;t) Fi{xi + l;t) 



E 

M-l<X2<xi 



E E 

M—l<X2<xi yi=xi 



Fo(x2 + l;t) Fo(yi + l;t) 



If we denote the determinant on the RHS by f{x2,yi), this is equal to 

OO OO J/l 



Since 



we have 



<X OO CO 



E E E E E E n-^'y^) 

X2=M-1 Xi=X2 + l yi=Xi X2=M~l yi=X2 + l Xi=X2 + l 



oo oo 



CO CO 



X2=M yi=X2+l X2=M xi=X2+l 



Fo(a;2 + l;t) F^^x^ + l-t) 

Fo{x2] t) - Fo{x2 + 1; t) Fo(xi; t) - Fo{xi + 1; t) 

^xi+X2+l ^-2t 



(x2+iy.'^ (xi+iy.^ 

Xll Xl\ 



OO CO 



(a;i + l)!(x2 + l)! 



,Xi -X2 



[x2{t)>M-l]= J2 Y i^2-x,) 



t 



XI+X2 + 1 



-2t 



X2=M+1 x\=X2+l 
oo oo 



^E E(---.) 



4^Xl+X2 

Xi\X2\ 



xi=M X2=M 

This is nothing but the expression (17. lip . (17. 12p for = 2. 



(7.13) 



(7.14) 



(7.15) 



(7.16) 



7.3 Position of particles: Step case 

In the previous subsection we considered the fluctuation of a single particle. To study the 
joint distribution, the original form of the Green's function is not very convenient. We 
rewrite it in a form such that the position of particles can be regarded as a dynamics of 
the first walker in a certain vicious walk [57,58]. 
Let us consider a determinantal weight. 



N-l 

n 

r=l 



detfC+iK, ^r%,i=i,...,r+i ■ det[^f ^(xf )],=o,...,^_^i, 



1=1,. ..,N 



(7.17) 
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under the condition < x\,x\ < x^ < X3, ■ ■ ■ and the convention ^^^.^^^(x^.,.]^; x^^^) = 1 
for r = 1, . . . , — 1, j = 1, . . . , r + 1. Here the functions ip'^p {x) and 0^ i^i u) defined 
by 



2Tii 



dw 



w 



r-l-j 



wt 



ro,- 



(1 + 



and 



In particular 



0^i,r-2(a;i,a;2) 



0^,.+l(xi,X2) 



Xi — X2 — 1 
r2 - ri - 1 



, Xi > X2, 

0, xi < X2. 



(7.18) 



(7.19) 



(7.20) 



This can be considered as a weight of a vicious walk(Fig. 10). A different feature from 
(15. 6p and (16. 7p is that the number of particles is not fixed in (17.171) . 
Let Pg denote the corresponding measure. Then we have 



G{xi, 
For N = 2 this reads 

G'(xi,X2;t) 



xn; t) = Pg[xi = xi,xl = X2, . . . ,x^ = xn]- 



(7.21) 



Fo{x2-y2;t) Fi{xi-y2]t) 
F^i{x2-yut) Fo{xi-yut) 



E 



1 1 



4^)(x?) 4'\xl) 



where and can be proved easily by using properties the determinant and 

F„ {x; t). The proof for general is similar and is given in [58]. The merit of rewriting this 
way is that the weight (I7.17P has a similar form to (I4.27P for tGUE. As a consequence we 
can describe the joint distribution by a Fredholm determinant. The kernel takes the form, 



KN2{ri,Xi;r2,X2) = K^2iri, Xi;r2, X2) - ^^^^^^(xi, Xa) 



G 



with 



N-l 



K^2iri,xi;r2,X2) = ^iljl'''\xi)(f^p\x2) 

j=0 



(7.22) 



(7.23) 



(r) (r) 

where (pj is a polynomial orthonormal to ip] . For the step initial condition for which 

yj+i = -h {3 = 0, 1, 



N — 1) can be taken as 



Vf\x) 



'l + z 



-.-zt 



To 



(7.24) 
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where Tq is the contour enclosing the origin anticlockwise. As a result we have a double 
integral expression of the kernel, 

J-P^ 2m 2m (1 + w)'^^+'^^{—zy^[w — z) 

By doing the asymptotic analysis as usual, we can reproduce the result in subsection 6.5 
for the cTASEP. 

7.4 Position of particles: Alternating case 

In the previous two subsections we have seen that the method using the Green's function 
allows us to reproduce the results for the step initial condition, which were obtained by a 
combinatorial approach. A big advantage of the method of using the Green's function is 
that one can deal with arbitrary initial configuration of particles at least up to some stage 
of the analysis. 

In particular in [57,58] the joint distribution for the alternating case was analyzed. For 
a technical reason we take the initial condition to be Uj+i = —2j {j = 0,1, . . . , N — 1). 
But if we take t finite and large enough, deep inside the negative region x{< 0) the 
dynamics of particles is the same as that for the alternating initial condition. Remember 
here that up to now the multi-layer PNG techniques have not been successfully applied to 
obtain the multi-point statistics for this initial condition. The discussions are the same as 
for the step case until (17.231) . The difference is that this time we employ a different set of 

(r) (r) 

(fi) s which are orthonormal to ip) s. The functions are found to be 



2m Jro ^ + z z' ^ 



3 



Then the kernel becomes 

K%^{ri,Xi]r2,X2) = K^i(ri, Xi; r2, X2) - 4>f^^,.^{xi,X2), 

f dw f dz {1 + zY^+'-^-^{-wY'{l + 2z)e^'"-'^^ 
K^,{n, X,- T2, X2) - ^ — J^^ — ^ ^)-i+n-i(_;,)..(^ _ + ^ + ^) • ^^-27) 

By applying the saddle point analysis, one gets a limiting kernel which is different from 
(I4.39p . Though (17.271) looks similar to (17.251) . the pole at z = —1 — w makes the situation 
completely different. 



8 Concluding Remarks 

In this paper we have discussed the applications of random matrix techniques to the study of 
fluctuations of the one- dimensional totally asymmetric simple exclusion process(TASEP). 
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Besides a brief review of the techniques, we have explained three approaches. The first one 
is the combinatorial approach, the second is the discrete polynuclear growth model and the 
last one is the Green's function method. The point is that in all these approaches there 
appear some objects whose weight can be described by a product of determinants and that 
the same kind of weight appears in the random matrix theory. We have mainly explained 
the methods and results for the step initial condition. Other cases including the alternating 
initial condition are also briefly explained. We have also seen the relationship of the TASEP 
to the directed polymer in random medium and the two types of the polynuclear growth 
models. 

In spite of its many successful apphcations, it has not yet been clear to what extent 
this method could be extended to study other properties of the ASEP. In fact there are a 
lot of open questions. For instance the correlation of the current at the origin N{ti), N[t2) 
at two times is a natural quantity to consider, but we don't know if it is handled by a 
generalization of the random matrix techniques at the moment. More studies with new 
ideas are certainly desired. 
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Appendix A 



In this appendix we explain some definitions, notations and formulas related to the Young 
tableaux. The partition A = (Ai, A2, . . .) is a collection of nonnegative integers Aj's which 
satisfy Ai > A2 > . . . and is conveniently represented by a Young diagram. For instance 



the partition A = (4, 3, 1, 1) corresponds the Young diagram 



A semistandard 



Young tableau(SSYT) is a Young diagram with each box filled with a nonnegative integer 
with the condition that they are non-decreasing in each row and are increasing in each 
column. We often restrict the entries of integers to fill to be from a finite set {1,2,..., N}. 
The Schur function sx{xi,X2, . . .) is defined by 



Sx{Xi,X2, . . .) 



T:SSYT 



where Tj is the number of i in SSYT. For instance when X — (2,1), X4 — x^ — . 
we take entries from {1, 2, 3}, all possible SSYTs are 



2 



J_ 
3 



2 



J_ 
3 



2 



J_ 
3 



_2_ 
3 



(A.l) 
. = and 

(A.2) 
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and hence 

sxixi,X2, xs, 0, . . .) = xlx2 + xjxs + Xix\ + Xixl + x\x^ + X2x\ + 2X1X2X3. (A.3) 
The number of the SSYT's of shape A with entries from {1,2,..., A^} is given by 

L(A,7V)=s,(l_^,0,...). (A.4) 

N 

There is another representation of the Schur function, 

det(x^"'+^-')i<,,,-<^ r 
det(xf"*)i<ij<iv 

By putting Xj = a^"^ {i — 1, . . . , N) and then taking the hmit a — > 1, one gets a formula. 



SA(a;i,...,Xiv,0, ...) = Tv-ix — ■ i^-^) 



..(i,...,i,o,...)= n ^-^i (A-6) 

' ' l<i<j<N ^ ^ 
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Figure Captions 

Fig. 1: The one-dimensional asymmetric simple exclusion process with parallel update. 
Fig. 2: The step initial condition for the TASEP. 

Fig. 3: An example of space-time trajectories of particles. The trajectories of the rightmost 
four particles are shown until each particle has made the fourth hops. 

Fig. 4: The dynamics of the eigenvalues of the GUE Dyson's Brownian motion (tGUE). 

Fig. 5: The waiting time table. The time coordinate s and space coordinate x in the 
discrete PNG model are also shown. 

Fig. 6: The discrete PNG model. The second layer in the multi- layer PNG model is also 
shown. Initially at time s = the surface is flat and nothing happens until s = 1. At 
time s = 1 a nucleation of height one occurs at the origin. Between time s = 1 and s = 2 
the surface grows laterally in both directions with unit speed. The height at the origin at 
s = 3~ is determined by the "higher one wins" rule. The height difference 1 becomes the 
height of a nucleation at the second layer. 

Fig. 7: A snapshot of a simulation of the multi-layer discrete PNG model. Note a similarity 
to the Fig. 4. 

Fig. 8: The alternating initial condition for the TASEP. 

Fig. 9: Another discrete PNG model in subsection 6.5. The height of the top layer gives 
the positions of particles through (16. 4p . Nucleations at a lower layer occur stochastically. 

Fig. 10: Vicious walk related to the Green's function of TASEP. The height of the top 
layer gives the positions of particles. The dependence on the initial condition is reflected 
in the function ipj{x). 
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